Level density of ^°Fe and low-energy enhancement of 7-strength function 
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The ^^Mn(d, n)^^Fe differential cross section is measured at Ed = 7 MeV. The ^''Fe level density 
obtained from neutron evaporation spectra is compared to the level density extracted from the 
^^Fe(^He,Q7)^^Fe reaction by the Oslo-type technique. Good agreement is found between the level 
densities determined by the two methods. With the level density function obtained from the neutron 
evaporation spectra, the ^®Fe 7-strength function is also determined from the first-generation 7 
matrix of the Oslo experiment. The good agreement between the past and present results for the 
7-strength function supports the validity of both methods and is consistent with the low-energy 
enhancement of the 7 strength below ~ 4 MeV first discovered by the Oslo method in iron and 
molybdenum isotopes. 



I. INTRODUCTION 

The unusual low-energy {E.y ^ 4 MeV) enhancement 
of the 7-strength function has been found recently - first 
for ^^Fe and ^^Fe nuclei [l|, Q and then for the set of 
molybdenum isotopes Q- The ('^He,a7) and (•^He,"^He'7) 
reactions and sequential extraction procedure developed 
at the Oslo Cyclotron Laboratory were used for this pur- 
pose. These results contradict the existing understanding 
based on different extrapolations of the giant dipole res- 
onance (GDR) to the low-energy region. The strength of 
the observed enhancement may indicate a different mech- 
anism (in contrast to the GDR) in the low-energy region 
of the 7-strength function. This requires additional the- 
oretical efforts to explain the observed enhancement. 

The Oslo method allows one to extract both the nu- 
clear level density (NLD) and the 7-strength functions 
from the first-generation 7 matrix P{Ei,E^) obtained 
from particle-7 coincidences in the (■^He,a7) reaction. 
Although it has been established that the method works 
reasonably well in practice, the question of the applica- 
bility of the Axel-Brink hypothesis remains open. This 
hypothesis assumes that the 7-strength function depends 
only on the energy of the 7 transition and not on the 
excitation energies of the initial Ei and final [Ei — E-^) 



states U, |5j] . This assumption leads to the factorization 
of the first-generation 7 matrix obtained from an Oslo 
experiment as: 



P{E,,E^) ^ p{E, ~ E^)T{E^) 



7/' 



(1) 



where p is the NLD and T is the radiative trans- 
mission coefficient, which is connected to the 7- 
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strength function through the relation: T{E^) — 

"^T^^XL fxL{E^)E\ . The p and T functions are 

determined by an iterative procedure j^ through the ad- 
justment of these two functions at each data point until a 
global y^ minimum with the experimental P{Ei, E^) ma- 
trix is reached. Another assumption is that the 7 tran- 
sitions originating from some energy interval feed levels 
with the same decay properties as those populated in 
the ('^He, a) reaction at the same excitation energy in 
the residual nucleus. This assumption has been partially 
supported by comparison of the results from two different 
reactions, namely ('^He,a) and (^He,'^He') populating 
the same residual nucleus. Although the Oslo method 
has been thoroughly tested, concern remains particularly 
about the validity of the Axel-Brink hypothesis. For ex- 
ample, the theory developed in Ref. p\ claims that the 7 
strength for spherical nuclei should depend on the tem- 
perature of the final state, implying that the Axel-Brink 
hypothesis is not valid. 

In order to address the above concerns, the NLD in 
Eq. (Q should be measured independently by a differ- 
ent kind of experiment. One of the most reliable meth- 
ods to extract the NLD below the particle separation 
threshold is based on measurement of particle evapora- 
tion spectra from nuclear reactions. Such spectra are 
described by a simple model of nuclear reactions based 
on Hauser-Feshbach formalism; according to this formal- 
ism the shape of the particle spectra depends only on 
the NLD of the final nuclei and the transmission coeffi- 
cients of outgoing particles. Because transmission coeffi- 
cients can be tested experimentally through the capture 
cross section of an inverse reaction, the NLD can be de- 
duced from spectra. The concern with this method is 
with the possible contribution of prc-cquilibrium and/or 
direct reaction mechanisms which potentially can distort 
the shape of statistical evaporation spectra and thus af- 



feet the NLD functions. The Oslo method uses 7 transi- 
tions which have been proven to be statistical. Therefore 
the comparison of NLDs obtained from these two ex- 
periments employing different reactions and extraction 
procedures will not only enable the comparison of the 
two methods, but also to estimate the possible distortion 
of particle evaporation spectra caused by direct reaction 
contributions. 

In this work the ^^Mn(d, n)^^Fe reaction was inves- 
tigated. The NLD obtained from the neutron evap- 
oration spectrum was analyzed and compared to the 
NLD determined from the Oslo experiment, where the 
^^Fe('^He,a7)^^Fe reaction was used. Since the method 
and results of the latter experiment have been thoroughly 
described in recent publications [^ |9j , we concentrate on 
the description of the ^^Mn{d,n)^^Fe experiment. 
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II. EXPERIMENT AND METHOD 

The experiment was performed with a 7 MeV deuteron 
beam from the John Edwards Accelerator Laboratory 
tandem at Ohio University. To measure the neutron spec- 
trum, the beam swinger facility 10] has been used. This 
allows the measurement of angular distributions by ro- 
tating the incoming beam and the target chamber with 
respect to the direction of outgoing neutrons. A self 
supporting 0.74 mg manganese foil was used as a tar- 
get. The energy of the outgoing neutrons is determined 
by the time-of-flight method with a 7 m flight path and 
NE213 neutron detectors. A 3 ns pulse width provided 
an energy resolution of about 100 keV and 800 keV for 
1 and 14 MeV neutrons, respectively. The neutron de- 
tector efficiency was determined with the calibrated neu- 
tron flux from the ^^Al(d, n) reaction on a thick Al tar- 
get at Ed = 7.44 MeV [ij- This allowed the determi- 
nation of the detector efficiency from 0.2 to 14.5 MeV 
neutrons with an accuracy of ^ 6%. The neutron spec- 
tra were measured at nine different angles from 20 to 150 
degrees to determine the angular distribution of outgo- 
ing neutrons. Additional measurements with an empty 
target were performed at each angle in order to deter- 
mine the background contribution. The absolute cross 
section was determined by taking into account the tar- 
get thickness, the accumulated charge of the incoming 
deuterons and the detector efficiency. The angular dis- 
tribution of outgoing neutrons is shown in Fig. ^ The 
observed angular anisotropy indicates the contribution 
from non-compound reactions at angles less than ~ 70°. 
The contribution of the non-isotropic part to the total 
neutron cross section is estimated to be about 30%. The 
deuteron break-up mechanism is also responsible for the 
cross section anisotropy at lower neutron energies. The 
angular dependence of the cross section at backward an- 
gles is flat and assumed to be due to the compound nu- 
clear mechanism. Therefore the spectra averaged over 
backward angles have been used to extract the NLD of 
the residual ^^Fe nucleus (see Fig. 12)). 



FIG. 1: The angular distribution for different energy groups 
of outgoing neutrons from ''''Mn(d, n) reactions. 
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FIG. 2: The neutron evaporation spectrum at backward an- 
gles (points). The spectrum calculated with the Empire pro- 
gram taking into account the level density according to the 
microscopic model of Ref. Llq| (solid line). The vertical line 
shows the neutron energy limit beyond which the first stage 
contribution to the total yield is greater than 90%. 



The procedure used to extract the NLD from the evap- 
oration spectra was proposed in Ref. 12]. This extrac- 
tion procedure is based on the Hauser-Feshbach theory 
of compound nuclear reactions, according to which the 
particle emission cross section is: 



3— (£a,eh) 

deb 



(2) 



ECN/ N E/TT ^b{U, J, TT, E, /, 7r)pb{E, I, tt) 



Jtt 



with 



r{U,J,7:)^J2lj2Th'iU,J,7r,Ek,Ik,7^k) + 

b' \ k 
nU-B,, 

Y, I dE'Tb'iU,J,Tr,E',I',Tr')pb>{E',I',7T') 



(3) 



Here g'~^^ {sa) is the fusion cross section, Sa and e^, are en- 
ergies of relative motion for incoming and outgoing chan- 
nels (eh = U — Ek — Bb, where Bi, is the separation energy 
of particle b from the compound nucleus) , the Tb are the 
transmission coefficients of the outgoing particle, and the 
quantities (f/, J, tt) and {E, I, it) are the energy, angular 
momentum, and parity of the compound and residual 
nuclei, respectively, Ec is the continuum edge. It follows 
from Eq. that the shape of the NLD is determined by 
both the transmission coefficients of outgoing particles 
and the NLD of the residual nucleus Pb{E, I, tt). Trans- 
mission coefficients can be calculated from optical model 
potentials usually based on experimental data for elastic 
scattering and total cross sections in the corresponding 
outgoing channel. This leaves the NLD as the only un- 
known parameter, which can be extracted from this equa- 
tion by using the experimental differential cross section. 
Details and assumptions of this procedure are described 
inRefs. [Hill. 

Neutron transmission coefficients are calculated from 
the optical model potentials taken from the RIPL-2 data 
base ji^j. Ten potentials have been tested. These are 
potentials based on global systematics such as given by 
Wilmore and Hodgson [13 as well as those obtained for 
the local mass range near A = 56 nuclei. All of them have 
been found to give the same result (the same shape of 
neutron evaporation spectra) within ~10% for 1-15 MeV. 
Finally, the potential of Wilmore and Hodgson has been 
adopted and 10% errors are added to the neutron trans- 
mission coefficients. In order to extract the NLD of the 
residual nucleus the following procedure was adopted: (i) 
The NLD model is chosen to calculate the differential 
cross section of Eq. (PJ. The parameters of the model 
were adjusted to reproduce the experimental spectrum 
as closely as possible, (ii) The input NLD was improved 
by binwise renormalization according to the expression: 



Pb{E, I, tt) = pb{E, I, 7r)input 



{da /dsb) mens 

(dcr/de f,)calc 



(4) 



The procedure described above is only correct when the 
main contribution to the differential cross section comes 
from the first stage of the nuclear reaction populating the 
residual nucleus of interest. In our case the second stage 
contaminations open up above particle separation ener- 
gies and come mainly from (d, an) , (d , pn) and {d,nn) 
reactions. But as proposed in Ref. [Ig, as long as these 
contributions are less than the total error of the extracted 
NLD, the energy interval chosen for the extraction of 
the NLD can be extended beyond the particle separa- 
tion threshold. Assuming 10% experimental errors (or 



more), the excitation energy interval where the second 
stage contributions do not exceed 10% is 0-11.5 MeV for 
the ^^Mn{d,xn) reaction (see Fig.|2J. 



III. LEVEL DENSITY OF ^''Fe 

The extracted NLD for the ^^Fe nucleus is shown in 
Fig. |3| along with the density of discrete low-lying levels 
(upper panel) as well as with the NLD obtained from the 
Oslo experiment (lower panel). The good counting statis- 
tics in the region corresponding to the location of known 
discrete levels allowed an absolute normalization of the 
extracted NLD. This is necessary because the scaling fac- 
tor also depends on NLDs and transmission coefficients 
of other outgoing channels (mainly proton and a). 

FigureOdemonstrates that up to about 6 MeV excita- 
tion energy the extracted NLD almost perfectly follows 
the shape of the NLD function based on discrete levels. 
This implies that the transmission coefficients used in 
calculating the theoretical spectrum are correct. Above 
this point the discrete levels are not complete and their 
density drops, while the NLD obtained from our experi- 
ment continues to increase. It also agrees well with the 
NLD obtained from the Oslo experiment l9l| . One can see 
the good general agreement up to ^ 8 MeV. 

Both curves show the same step structures at around 4 
and probably at 6 MeV excitation energy. Similar steps 
have been interpreted in Q as a result of the breaking 
of Cooper pairs. Above about ^ 7.8 MeV these curves 
start to diverge, ending with differences of about 50% 
at an excitation energy of ~ 8.8 MeV. This deviation 
may stem from imperfections of the two methods, per- 
haps connected with violations of the underlying assump- 
tions. Some local deviations of the neutron transmission 
coefficients are also possible. However, the almost per- 
fect agreement below ^ 8 MeV gives confidence in both 
methods. 

Generally, one can conclude that in spite of the fact 
that these two methods use different underlying assump- 
tions, different nuclear reactions, and different mathe- 
matical techniques to extract NLD, a fairly consistent 
result has been obtained. This implies that the statis- 
tical mechanism dominates in both reactions. The ob- 
served difference between NLDs at higher excitation en- 
ergy requires further investigation. The observed small 
bump in the differential cross section at '^ 9 MeV of the 
neutron energy transforms to the corresponding bump in 
the extracted NLD at about 5.5 MeV of excitation en- 
ergy. Previously, similar structures were observed in Oslo 
type experiments for a variety of nuclei and have been ex- 
plained as a result of pairing correlations. The NLD of 
^^Fe from the Oslo experiment also exhibits structure at 
~ 5.5 MeV, although the shape is slightly different from 
that observed in the particle evaporation spectrum in Fig. 
O This is probably due to systematic uncertainties still 
inherent in both methods. 

The presence of the step structure in the NLD is an im- 
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FIG. 3: The comparison of the NLD extracted from neutron 
evaporation spectra (full circles) with discrete NLD (upper 
panel) and with NLD (open circles) obtained from Oslo-type 
experiment (lower panel). 
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FIG. 4: The NLD extracted from neutron evaporation spectra 
(full circles) compared to Fermi-gas (dashed) and microscop- 
ical model (open circles) calculations. The full line shows the 
fit to the data. 



portant feature from a practical point of view, because 
this might introduce corrections to available systematics 
of NLD parameters widely used for calculation of reac- 
tion rates. Almost all systematics are primarily based on 
NLDs obtained from neutron resonance spacings at the 



neutron separation energy J5„ and the NLD of low-lying 
discrete levels. To calculate the total NLD at Bn the spin 
cutoff parameter is used; this parameter is not known to 
high accuracy in this energy region. The NLD in the in- 
termediate region is often obtained by a simple interpo- 
lation between these two anchor points, assuming some 
NLD model. The Fermi-gas model is most often used. 
For example the NLD according to the Fermi-gas model 
with parameters from the latest systematics of Ref. [13 
is shown in Fig. ^ This curve fits the discrete NLD well, 
but overestimates the experimental points by a factor 
of ~ 1.7 at higher excitation energies. A slightly bet- 
ter result is given by the microscopic model ,18.] recom- 
mended in the RIPL-2 data base 0|. However, because 
this NLD is also renormalized to the neutron resonance 
spacing and density of low-lying levels, the discrepancy 
is still present. These or other similar models may cause 
a sizable overestimation of calculated reaction cross sec- 
tions and incorrect ratios of cross sections for different 
channels. 

To fit the experimental NLD the Fermi-gas formula has 
been adopted: 



piE) 



exp{2y/a{E - S)) 



(5) 



12V2ai/4(^_^)5/4^' 

where the spin cutoff parameter a based on a rigid body 
moment of inertia is expressed as: 



a^ = 0.0U5A^^^^/{E-S)/a. 



(6) 



The Fermi-gas parameters which give the best fit to 
experimental points between 2 and 11 MeV are a — 
5.65(10) MeV-i and S = 0.65(20) MeV. The spin cut- 
off formula based on half of the moment of inertia does 
not cause large changes. These parameters are slightly 
smaller than those obtained from the systematics of Ref. 
[13 (a = 6.16 MeV-\ S = 0.8 MeV). However, this 
difference is enough to cause ~ 70% discrepancy in cor- 
responding NLDs at 6-10 MeV excitation energies (see 
Fig.EJ). 



TABLE I: NLD parameters of ^®Fe from different systematics. 



Systematics 


a 


S 


T. von Egidy 


6.16 


0.8 


G. Rohr 


5.61 


2.81 


S.I. Al-Quraishi 


5.98 


1.68 


Present 


5.65(10) 


0.65(20) 



The NLD parameters for ^^Fe have also been obtained 
in Ref. il9l| and uu\ from {p,p'), {p,a) and {a, a') reac- 
tions. Despite the fact that the same reactions and the 
same techniques have been used in these experiments, dif- 
ferent results were obtained. The first measurement re- 
ported a = 5.7 MeV~^ and 5 — 0.7 MeV (in good agree- 
ment with our results), while the second measurement 
gives a = 6.5 MeV~^ from {p,p') and a = 7.0 MeV~^ 



from {p, a) reactions. Part of the problem might be 
the strong correlation of the a and i5 parameters when 
they are extracted from particle evaporation spectra. In 
the case when the extracted NLD is not normalized, one 
needs to assume some value of S in order to obtain the 
value of a. As has been shown in [l3|, a 1 MeV increase of 
6 causes approximately 1 MeV~^ decrease of parameter 
a (the correlation is negative). 

More recently V. Mishra et al. ^^ reexamined this 
problem. Their result confirmed the conclusion of [l9j 
if it is assumed that the slope of the NLD is matched. 
On the other hand, if the parameters are extracted from 
absolute value of NLD at a particular point (as in neutron 
resonance analysis), the correlation between a and S has 
positive sign. 

In addition to the von Egidy compilation, other compi- 
lations have been prepared by Rohr 22] and Al-Quraishi 
et al. |23|- These parameters are shown in Table I. It 
is interesting to note that the three compilations have 
predicted slopes which are within 10% of one another 
(same nuclear temperature). The magnitude predicted 
for the NLD between 6 and 10 MeV differs by more than 
an order of magnitude. Rohr predicts the smallest NLD, 
Al-Quraishi next highest, the present data are second 
largest and von Egidy et al. predict the highest NLD. 
These differences in magnitude point out the importance 
of performing absolute magnitude normalization in in- 
ferring NLD. In this work, we were able to obtain the 
absolute NLD function by normalizing it to the NLD of 
known low- lying levels. This gives us the possibility of 
obtaining both the a and 6 values simultaneously. 

Thus, in order to establish the systematics of NLD pa- 
rameters, it is necessary to take into account the absolute 
NLD for the whole excitation energy interval. Such infor- 
mation can be obtained either from particle evaporation 
spectra by using Eqs. (PJ and Q or from the Oslo type 
experiments by using the sequential extraction method 
applied to the particle-7 coincidence matrix. 



IV. 7-STRENGTH FUNCTION IN ^'^Fe 

The 7-strength functions for ^^Fe and ^^Fe have been 
obtained from an Oslo type experiment in Ref. [3|. The 
striking feature of these functions is the increase of 7 
strength in the region below 4 MeV; this cannot be under- 
stood within existing models. The main concern about 
the Oslo method is that possible violation of the Axel- 
Brink hypothesis might result in some systematic devia- 
tion of the obtained 7 strength - that is, if the 7-strength 
function depends not only on the energy of the 7 transi- 
tion, but also on the temperature of the final state. Such 
temperature dependence may stem from the temperature 
dependence of the GDR width caused by different damp- 
ing mechanisms debated in the literature [2J| . According 
to the Fermi-liquid model |2ItI | the width is determined by 
the collision of quasi particles in the nuclear volume that 



results in the following temperature dependence: 

TiE^,T)^iE?^+ATr^T^), (7) 



where T — \JU jo). But as has been shown in Ref. |2q . 
the systematic deviation of the 7-strength function due 
to such temperature dependence is about only 15% in the 
region below 2 MeV. The other mechanism is connected 
to nuclear shape fluctuations, leading to a square root 
temperature dependence of the GDR width. The exper- 
iments on Sn and Pb [23 support the temperature de- 
pendence of the GDR, but the mechanism responsible for 
such effects is still under debate. A good fit to these ex- 
perimental data is obtained with the Fermi-liquid model 
accounting for damping of the GDR according to Eq. |Zj), 
while taking into account the dipole-quadrupole interac- 
tion term arising from the nuclear deformation |23| . 

The low-energy enhancement of the 7-strength func- 
tion observed in the Oslo experiment for iron isotopes 
might stem from different (i.e. not GDR) modes of nu- 
clear excitation and therefore might exhibit a different 
temperature dependence. Therefore it is important to 
investigate the validity of the Axel-Brink hypothesis in 
the Oslo method for the iron isotopes. For this purpose 
we can use the first generation matrix P(£'i, E^^ obtained 
from the Oslo ^^Fe('^He,a7)^^Fe experiment and the NLD 
from ^^Mn((i, n)^^Fe reaction to obtain the 7-strength 
function of the ^^Fe, giving 



f{Ey,Ei) — 



1 N{E,)P{E,,E^) 
2tt p{E,~E-,)E^ ' 



(8) 



where Ei = E^ + Ef. It is clear that the obtained 7- 
strength function should not deviate considerably from 
that extracted solely from the P{Ej, Ej.) obtained by the 
sequential iterative procedure because the corresponding 
NLD functions agree well (see Fig. |21l. However, due to 
slightly different slopes of these functions in the region 
of 4-7 MeV and ~ 1.7 times disagreement above this re- 
gion (which in principle can be caused by temperature 
effects), it is interesting to investigate the magnitude of 
corresponding local changes in the 7-strength function. 
The normalization constant N{Ei) in Eq. © has been 
determined in Ref. [3 at i^i = Bn by using supplemen- 
tary experimental information from systematics of total 
7 widths of initial states (neutron resonances). The com- 
parison of 7-strength functions obtained from both the 
sequential extraction Oslo method and Eq. ^ using the 
NLD from the evaporation spectrum is shown in Fig. 
|31 There is no significant disagreement. It can be con- 
cluded that possible temperature effects on the extracted 
7-strength function are rather small compared to total 
uncertainties in the experimental data. Thus the appli- 
cability of the Axel-Brink hypothesis in the Oslo method 
is justified within the accuracy of the experimental data. 
However, one should keep in mind that in both cases 
investigated here, one extracts an averaged 7-strength 
function for a wide region (several MeV) of final energies. 
In order to detect a possible temperature dependence of 
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FIG. 5: The 7-strength function obtained from Oslo first- 
generation matrix P{E^,Ex) with NLD from (d, n) reaction 
(filled circles). The 7-strength function obtained solely from 
P{Ej,Ex) (open circles). 



the 7-strength function, one could extract the 7-strength 
function from several limited regions of final energy, how- 
ever, for this, one needs to increase the statistics of the 
experiment to be able to reduce statistical errors. 



a broader range of nuclei. This agreement helps to elim- 
inate most of the potential systematic errors inherent in 
these methods, including such important problems as the 
unknown contribution of direct processes in particle evap- 
oration spectra. The neutron evaporation spectrum from 
the (d, n) reaction measured at backward angles does not 
contain such a contribution. 

The NLD function of ^^Fe can be fit by the conven- 
tional Fermi-gas model in the region of 2-11 MeV exci- 
tation energy. However, a local deviation of the order of 
^ 40% has been observed at 5 MeV. The presence of this 
structure leads to disagreement in the Fermi-gas param- 
eters obtained from our current data and from available 
systematics. In order to verify the character of these 
structures, more experimental data are needed for neigh- 
boring nuclei. 

The 7-strength function for the ^^Fe isotope obtained 
in Ref. [J has been extracted by using the NLD from 
the neutron evaporation spectrum. The new 7-strength 
function agrees well with the previous one within exper- 
imental errors. This indicates small temperature effects 
on the extracted 7-strength function. Thus the Axel- 
Brink hypothesis used in the Oslo method appears to be 
justified within the present experimental uncertainties. 
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V. SUMMARY AND CONCLUSIONS 

The ^^Fe NLD has been extracted from the neutron 
evaporation spectrum of the ^^Mn(d, n)^^Fe reaction. 
This NLD has been compared to that obtained from the 
^^Fe('^He,a7)^^Fe reaction by using the sequential extrac- 
tion technique developed at the Oslo Cyclotron Labora- 
tory. The NLDs obtained from these two different types 
of experiments are in good agreement with each other. 
This indicates the consistency of these two methods and 
possibility of applying such methods to investigation of 
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